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Abstract 

Using a scheme based on a Mach-Zehnder interferometer, we propose an analysis of the super- 
position of polarized laser beams at a given angle. The focus of our study is the spatially varying 
polarization state of the resulting field, also known as a polarization grating, generated by this 
setup. Our proposal combines a theoretical description of the resulting field in terms of its Stokes 
parameters with an experimental demonstration of the existence of such a polarization grating due 
to the effects of polarization on beam interference experiments. 
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I. INTRODUCTION 



Understanding interference has been seminal in optics. More than two centuries ago, 
Young presented his Bakerian Lecture which contained an experimental demonstration of 
the general law of interference of light^. Fifteen years later, Fresnel and Arago studied the 
effect of the polarization state of light beams in the phenomena of interference. Thus, 
through interference, evidence of the transverse wave nature of light was brought forward. 

In the second half of the twentieth century there were many studies on the interference 
of polarized light for the undergraduate laboratory. Various interferometric methods were 
proposed to carry out such experiments using (i) a Young double slit experiment, covering 
both slits with different polarizing filters^i^i^; (ii) the polarizing properties of the ordinary 
and extraordinary axis in a double refracting crystal of calcite^*^ or a nematic liquid crystal^; 
(iii) a Ronchi grating as beam splitter placing polarizing filters afterwards the grating and 
analyzing at the conjugate plane of the grating through a focusing element^; (iv) a Mach- 
Zehnder interferometer with polarizing filters at each arm^°. Just to list a few schemes. 

In the same period of time, a sound mathematical description for the Young double 
slit experiment with polarized light was presented characterizing the resulting field by the 
Stokes parametersii. Also, an alternate mathematical formulation and a proposal for using 
this scheme in image processing involving Wollaston prisms was discussed^^. Recently, this 
was used to understand the concept of a quantum eraser, where the polarization of two 
macroscopic fields is manipulated to loose or restore interference fringes^. 

In all of the aforementioned studies, discussions focused on understanding the resulting 
intensity pattern, which is directly related to the Fresnel-Arago laws. To our knowledge, 
little has been said about the polarization state of the resulting field; it is possible to obtain 
a ^-linearly polarized field from the superposition of two collinear fields with right and 
left circular polarization, with the given angle 6 G (—90, 90] degrees defined by the phase 
difference between the superposing fields^. 

We believe the analysis of the superposition of two light fields could go beyond the study 
of Fresnel-Arago laws. Our motivation comes from the fact that a rich spatially dependent 
polarization structure, known as a polarization grating, arises from the superposition of two 
non-collinear polarized light beams with different polarizations. This polarization grating 
can be used to simplify the fabrication of three dimensional periodic microstructres elimi- 
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nating complicated procedures to optimize contrast and polarization^^. It can be recorded in 
polarization sensitive materials, like bacteriorhodopsin, to perform phase shifting shearing 
interferometry without neither mechanical moving parts nor mounted gratings^^. It can also 
be used to simultaneously rotate or control the orientation of multiple microscopic birefrin- 
gent particles, which might be useful for biomedical applications^^. These are just a few 
examples involving polarization gratings in modern research. 

Our experimental proposal uses a slightly modified Mach-Zehnder interferometer to pro- 
duce light with spatially varying polarization. Our theoretical description of the experiment 
is based on Jones calculus^ and Stokes parameters^ for polarized monochromatic light. 
The analysis is complemented with three specific examples involving the superposition of 
combinations of linearly and circularly polarized light showing some specific polarization 
gratings that can be produced with the proposed scheme. The presence of such polarization 
structures is experimentally confirmed through qualitative analysis of the resulting light field 
with a linear polarizer. We encourage the reader to reproduce the experiment and calculate 
the values for the Stokes parameters of the resulting field^"^. 

For readers who may not be familiar with Jones calculus and Stokes parameters we suggest 
the textbooks by Simmons and Guttmann^S or by Born and Wolf^^ in addition to the articles 
mentioned above^^^^^. For those interested in an advanced treatment, we recommend as 
a starting point the articles by Tervo et. af^ and Roychowdhury and Wolf^^, as well as the 
introductory book by Wolf2^. Their analyses deal with polarization and coherence degree of 
superposed arbitrary electromagnetic fields in three dimensional space. 

II. THEORETICAL ANALYSIS 

Figure [TJillustrates a simplified version of the superposition scheme. We focus on the plane 
of incidence defined by the xz-plane so the y-coordinate will be obviated. Two polarized 
monochromatic plane waves of light intersect with a small angle 6, such that sin 6' ^ 6' in 
radians, at some point p{x, z) = xx + zz on the detection line S. Such light fields are 
described by the equations 

E2{x,z,t) = ^2 e<'^'^^("'^)-"*+^^)i?,(^)£2(«2,52). (1) 



3 



where the distances di{x, z) are the distances from the i-th beam source to the point p{x, z), 
e.g. di{x, z) — z. The counterclockwise rotation of the polarization state of the second beam 
about the y-axis is introduced by means of the rotation matrix 



^ cos^ sin^^ 



Ry{0) 



1 
sin^ cos^ 



(2) 



The unitary polarization vector state ij{oLj, Sj) is, up to a phase constant, a Jones vector 

ij{aj, Sj) = cosctj X + e*'^-' sinctj y, (3) 

with parameters in the ranges aj G [0,7r/2] and Sj G (— tt, tt]. The symbols x and y are the 
unitary vectors in the x- and y-directions. 

The Stokes parameters for the total field E{x, z, t) — Ei{x, z, t) + E2{x, z, t), at a point 
p{x, z) on the detection line E are given by the expression 



Si = {E{x,z,t),aiE{x,z,t)) 

= sl'^El + sf^El + 2ErE2 Re [e'^*el • aiRy{e)i2] . 



(4) 



The angle brackets are shorthand notation for time averaging over the detection interval, 
which is large compared to the period associated with the optical radiation frequency, 
(w(r, i), v(r, i)) = |i*(r)* • v{r) for plane waves (asterisk meaning complex conjugation). 
The symbol Oi for i = 0, 1, 2, 3 denotes the Pauh matrices 
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The parameters s'f'' are the Stokes parameters for the j-th polarization vector ij{aj,6j) 

= i* ■ (Toij = 1, 



3 "u-j 

= i* • aiij = cos2aj, 

— i* • (72£j = sin 2a j cos Sj, 

— i* ■ a^ij ~ sin 2a j sin 8j. 



(6) 
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Finally, the phase difference parameter A$ can be approximated as 



= k X sin 6 + A0 

k X 9 + A(f), (7) 

with the initial phase difference between the sources of the beams given by A0 = 02 — 0i- 
In our experimental scheme the source for both beams is the same laser so the initial phase 
difference is null, A0 = 0. The experimental setup uses small angles, 6 < 10^^ radians, 
that for practical purposes, cos 6 ^ 1 — 6"^ and sin^ ^ 9. These experimentally feasible 
restrictions allow us to consider the polarization state of the second field beam E2 in its 
own reference frame and in the general reference frame almost equal; i. e. the Jones vectors 
^2(^2, (^2) and Ry(9)i2{oi2, 82) are almost parallel, 

(£2(«2, ^2), -R{/(^)^2(«25 (^2)) = cos 6' + ? sin ^ sin 2 a sin (5, (8) 

^ 1 -^^ + «^sin2asin(5, ^ < 10"^ 
~ 1. 

These approximations simplify the theoretical treatment, thus the real parts involved in the 
last term of Eq. (jl]) are given by 

Re (e*^*£^ ■ (To-Rj,(^^)s2) ~ cosai cosa2 cos A$ + sinai sin 02 cos (A$ + A(5) , 

Re (e*^*£* ■ aiRy{6)i2) ~ cos ai cos 02 cos A$ — sin ai sin 02 cos (A$ + A6) , 

Re (e*^*£* ■ a2Ry{0)i2) ~ cosai sin 02 cos (A$ + ^2) + sinai cosa2 cos (A$ — 5i) , 

Re {e"^^ £{ ■ a^Ry{9)i2) ~ cosai sin a2 sin (A$ + ^2) - sin ai cosa2 sin (A$ - . (9) 

As usual, the Stokes parameter Sq is useful for discussing the intensity profile at the detection 
line such as discussed by Pescetti^ or CoUetii, while the latter three parameters, 5*1 to 6*3, 
relate to the polarization state of the field. 

Our purpose is to understand the polarization properties of the total field. In order 
to do so, let us consider the interfering beams carrying orthogonal polarizations, that is 
(£i(q;i, £2(a2, ^2)) = 0. Notice that two orthogonal polarization vectors can be written 
as ii{a,5) and i2{ot — 7r/2,5) if we relax the restrictions on the domain of a. The Stokes 
parameters for orthogonal polarization vectors fulfill the condition sf'' = —s\^^ for i = 1,2, 3; 
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the points on the polarization sphere that represent these vectors being antipodes. It is also 
possible to parametrize the amplitudes of the fields as /3 = arctan in a range P E [0, 7r/2], 
such that the corresponding normalized Stokes parameters for the total electromagnetic field 
on the detection line are 

So ~ 1, 

S^i ~ cos 2a cos 2/? + sin 2a sin 2/5 cos A$, 

5*2 ~ sin 2a cos 2/5 cos (5 — sin 2/9 (cos 2a cos (5 cos A$ — sin (5 sin A$) , 

S3 ~ sin 2a cos 2/3 sin 5 — sin 2/3 (cos 2a sin 5 cos A$ + cos (5 sin A$) . (10) 

It is possible to write the latter three normalized Stokes parameters, 5*1 to 5*3, as the vector, 



S = Si si + S2 S2 + S3 S3 

^ i?,,(7r -5)i?,3(2a) (sin2/3 ^+ cos 2/3 Si) , 

with the vector g defining a great circle on the S2S3-plane of the polarization sphere. 



g = cos A$ S2 + sin A<l> Sg 

and the rotation matrices given in the traditional way 

/ 



(11) 



(12) 



1 

COS ^ sin -(9 
— sin i) cos i) 



RsM 



\ 



(13) 



cos'd sin'd 
— sin cos i} 

ly 

Equation [11] implies that the parameter a generates a counterclockwise rotation around 
the Sa-axis, the parameter (3 acts as a scaling factor and a si-translation on the great circle 
g, and the parameter 5 as a counterclockwise rotation around the Si-axis. 
The counterclockwise rotations were expected. The great circle g is obtained from the 
superposition of the fields emitted by two sources with equal amplitudes of emission and 
horizontal/vertical linear polarizations. The rotation -Rs^(7r — 5)-Rs3(2a) transforms the two 
sets of Stokes parameters that map the orthogonal polarization pair £i(0, 0), £2(^/2, vr) into 
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those mapping any other orthogonal pair ii{a,S) and £2(« — 7r/2,(5). Figure [2] shows an 
example of the effect of the set of parameters {a, (3,6} on the behavior of the polarization 
for the total field at the detection line. An interactive demonstration, where the user can 
input any combination of these parameters and obtain the corresponding Stokes parameters 
on the polarization sphere, is provided online^I. 

III. EXPERIMENTAL SETUP AND RESULTS 

We present the experimental realization and discussion of three cases that can shed more 
light on the problem when working in the undergrad laboratory. In the first two cases, beams 
are used with polarization states orthogonal to each other, Eq. (fTT|) . In the third case, the 
general treatment is used, Eq.(jl]). 

The experimental setup is shown in Fig. [31 Using a Mach-Zehnder interferometer, two 
beams are superposed at the back aperture of a microscope objective. The image of the 
beams superposition is formed at the focal region of the objective. Each of the beams is 
given a specific polarization state through a suitable retarder at the corresponding arm of 
the interferometer. Characterization of the superposition polarization state is performed by 
placing a linear polarizer as analyzer behind the focus of the objective. Images are captured 
for angles of 0, 7r/4, tt/2, and 3tt/A radians of the linear polarizer axis with respect to the 
horizontal axis. 

A continuous wave solid state laser, emitting at a 532 nm wavelength with a linear vertical 
output polarization state, ratio 100:1, is used as a source. The beam splitters, BSl and BS2, 
are non-polarizing and one of them, BS2, is mounted on a linear displacement stage in order 
to control the angle of interference, 6 in Fig. [H In the interferometer, each beam reflects 
once from a beam splitter and once from a mirror, thus canceling the phase shift introduced 
by each single reflection. A 40x microscope objective is used as an imaging element. A 
black and white l/2"-CCD, located just after the focal distance of the imaging element, is 
used to capture the images. 

The theoretical results shown here are calculated to fit the results at the CCD, whose 
range of detection is approximately x G [—6.5, 6.5] mm. An intersection angle given by 
6 ^ 6.1728 X 10~^ radians is used. 

A general interactive demonstration where the user can set at will all of the two beams 



7 



parameters is provided online^^. This demonstration generates the corresponding Stokes 
parameters on the detection hne, their representation on the polarization sphere, and the 
polarization ellipse for a given position on the detection line. 

In the following paragraphs polarization states will be denoted using the notation: cir- 
cularly right/left, TZ/ C, elliptically right/left, £TZ/£C, linearly 9, V{9). 

A. P-S Configuration. 

Balanced horizontal/ vertical linear polarization 

ai = 0, a2 = vr/2 , /? = 7r/4, 5i = 0, ^2 = t^- 

In this case, PI is a half- wave plate, whose fast axis is placed at 7r/4 radians with respect to 
the X-axis, producing horizontal linear polarization and P2 is removed to keep vertical linear 
polarization. Figure H] shows the intensities captured by the CCD camera at the four analyzer 
orientations for this configuration alongside the corresponding theoretical intensities. The 
vertical lines crossing both columns of figures are presented as reference marks to relate 
experimental with theoretical results. 

Figure [5] shows that the analysis is consistent with a polarization of the resulting field 
varying on a meridian of the polarization sphere crossing the ±45° linear polarization states; 
i.e. polarization varies periodically with the cycle: TZ P(45) £ — > P(— 45°) — » 7?., with 
intermediate states of elliptical polarization. 

The alphabetical labels (A,B,C,D) in both figures identify the intensities of the analyzed 
polarization states in Fig. HI as well as the related Stokes parameters, their corresponding 
mapping on the polarization sphere, and the polarization state they represent in Fig. [5l 
for the field at a given point x. Thus, the labels A/B/C/D show where the field has a TZ/ 
P(-45°)//:/P(45°) polarization state. 

B. R-L Configuration. 

Balanced right/left circular polarization 

ai = a2 = 7r/4 , j3 = 7r/4, 5i = —82 = vr/2. 

In this configuration, PI and P2 are quarter- wave plates whose fast axes are placed at 
±7r/4 radians with respect to the x-axis, thus their fast axes are perpendicular to each other. 
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producing right circular and left circular polarizations. Figure [6] shows the experimental and 
theoretical intensities corresponding to the four analyzer positions. The vertical lines are 
used to relate both results. 

Figure [7] shows that the analysis is consistent with a polarization of the resulting field 
varying on the equator of the polarization sphere; i.e. polarization is always linear with 
orientation angle varying periodically from —it to vr radians. It is shown that "The super- 
position of the right and left circularly polarized light yields linearly polarized light but the 
direction of the polarization depends on the phase angle between the two beams"—. 

The alphabetical labels in both figures identify the intensities of the analyzed polarization 
states in Fig. El as well as the related Stokes parameters, their corresponding mapping 
on the polarization sphere, and the polarization state they represent in Fig. [TJ for the 
field at a given point x. Thus, the label A/B/C/D shows where the field has a V{0°)/ 
P(45°)/P(90°)/P(-45°) polarization state. 

C. R-S Configuration 
Balanced horizontal linear and right circular polarization 

ai = 0, 02 = vr/4, (3 = 7r/4, 6i = 0, 82 = vr/2. 

Finally, here PI is a quarter wave plate, whose fast axis is placed at 7r/4 radians with 
respect to the a;-axis, producing right circular polarization and P2 is a half-wave plate, 
whose fast axis is placed at 7r/4 radians with respect to the x-axis, producing horizontal 
linear polarization. Figure [8] shows the intensities captured by the CCD camera at the four 
analyzer orientations alongside the corresponding theoretical intensities. 

Figure [9] shows that the analysis is consistent with a polarization of the field varying on 
some circle on the polarization sphere; i.e. polarization varies periodically being elliptically 
polarized but for two points where it is linearly ±45° polarized. 

The alphabetical labels show elliptical polarization states with varying eccentricity and 
inclination of the major axis, for this configuration. Their intensities. Fig. [8] and their 
related Stokes parameters, their corresponding mapping on the polarization sphere, and the 
polarization state they represent. Fig. [9], for the field at a given point x 
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IV. CONCLUSION 



We have presented an experimental scheme that an undergraduate student can use for 
analyzing the polarization state of the superposition of two slightly non-collinear polar- 
ized light beams. The equations modeling the Stokes parameters for this experiment have 
been presented. The explicit case of interfering orthogonal polarizations was discussed and 
complemented with two particular configurations to help elucidating this scheme; a third 
experimental configuration involving a general case, the interference of two non-orthogonal 
polarization beams, was also presented. It has been shown that the polarization state of 
light is spatially dependent in all cases due to the spatially dependent phase between the 
beams introduced by the impinging angle between them. 

Our experimental scheme can be implemented into an optical tweezer to demonstrate 
the transfer of intrinsic angular momentum to birefringent particles^- using polarization 
structures'^. This could also attract the attention to polarization, interference and mechan- 
ical properties of light at the undergraduate level. 
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Figure 2: (Color online) Effect of the parameters {a, /3, 6} for electromagnetic fields with or- 
thogonal polarization, A<I> G [0, 27r). (a) Linear polarization a € {0, vr/lO, 7r/5, Svr/lO, 27r/5}, 
/3 = tt/4:, 5 = 0; a = Oin dot dashed red. (b) Linear polarization a = 0, /3 G 
{7r/16,7r/8,37r/16,7r/4,57r/16,37r/8,77r/16}, 5 = 0; /3 = 7r/4 in dot dashed red.(c) Elliptical polar- 
ization a = 7r/4, /3 = 7r/4, (5 € {0, vr/S, 27r/5, Svr/S, 47r/5}; 5 = in dot dashed red. 
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Figure 3: Experimental setup. 
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Figure 4: (Color online) P-S Configuration. Interference of beams with horizontal and vertical 
linear polarization, equal field amplitudes and initial phases, ai = , 02 = '7r/2, 61 = 0, 62 = tt, 
P = 7r/4. First column presents the experimental intensities obtained after the analyzer. Second 
column present the theoretical intensities. The orientation of the analyzer corresponds to (a) 
horizontal, (b) 45°, (c) vertical, (d) —45°. The red vertical lines are presented as markers relating 
experimental and theoretical results. 
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Figure 5: (Color online) P-S Configuration. Interference of beams with horizontal and vertical 
linear polarization, equal field amplitudes and initial phases, ai = , 02 = vr/2, 61 = 0, 62 = vr, 
P = 7r/4. (Left) Stokes parameters Sq (solid black), Si (dashed blue), 5*2 (dot dashed red), 
5*3 (dotted green). (Top Right) Polarization trajectory on the polarization sphere given by the 
normalized Stokes parameters 5i, 6*2, 5*3. (Bottom Right) Polarization state corresponding to the 
point labeled A/B/C/D. 
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Figure 6: (Color online) R-L Configuration. Interference of beams with right circular and vertical 
linear polarization, equal field amplitudes and initial phases, ai = 02 = tt/A, P = tt/4, 61 = —62 = 
tt/2. First column presents the experimental intensities obtained after the analyzer. Second column 
present the theoretical intensities. The orientation of the analyzer corresponds to (a) horizontal, 
(b) 45°, (c) vertical, (d) —45°. The red vertical lines are presented as markers relating experimental 
and theoretical results. 
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Figure 7: (Color online) R-L Configuration. Interference of beams with right circular and vertical 
linear polarization, equal field amplitudes and initial phases, ai = 02 = 7r/4, /? = 7r/4, 5i = 
—82 = vr/2. (Left) Stokes parameters Sq (solid black), Si (dashed blue), 5*2 (dot dashed red), 
5*3 (dotted green). (Top Right) Polarization trajectory on the polarization sphere given by the 
normalized Stokes parameters Si, S2, S3. (Bottom Right) Polarization state corresponding to the 
point labeled A/B/C/D. 
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Figure 8: (Color online) R-S Configuration. Interference of beams with right and left circular 
polarizations, equal field amplitudes and initial phases, ai = , 02 = vr/4, 61 = 0, 82 = vr/2, 
El = E2- First column presents the experimental intensities obtained after the analyzer. Second 
column present the theoretical intensities. The orientation of the analyzer corresponds to (a) 
horizontal, (b) 45°, (c) vertical, (d) —45°. The red vertical lines are presented as markers relating 
experimental and theoretical results. 
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Figure 9: (Color online) R-S Configuration. Interference of beams with right and left circular 
polarizations, equal field amplitudes and initial phases, ai = , Q2 = vr/4, 5i = 0, 62 = vr/2. 
El = E2. (Left) Stokes parameters Sq (solid black). Si (dashed blue), 5*2 (dot dashed red), 
53 (dotted green). (Top Right) Polarization trajectory on the polarization sphere given by the 
normalized Stokes parameters Si, 52, 5*3. (Bottom Right) Polarization state corresponding to the 
point labeled A/B/C/D. 
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